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Abstract 

We study an integrable quantum field theory of a single stable particle with 
an infinite number of resonance states. The exact 5"-matrix of the model is 
expressed in terms of Jacobian elliptic functions which encode the resonance 
poles inherently. In the limit I — > 0, with I the modulus of the Jacobian elliptic 
function, it reduces to the Sinh-Gordon 5-matrix. We address the problem 
of computing the Form Factors of the model by studying their monodromy 
and recursive equations. These equations turn out to possess infinitely many 
solutions for any given number of external particles. This infinite spectrum of 
solutions may be related to the irrational nature of the underlying Conformal 
Field Theory reached in the ultraviolet limit. We also discuss an elliptic 
version of the thermal massive Ising model which is obtained by a particular 
value of the coupling constant. 



1 Introduction 



The search for solvable but realistic models in theoretical physics has always resulted 
in valuable by-products. In this respect, the subject of two-dimensional relativistic 
Quantum Field Theories (QFT) has played a prominent role in the last decade, since 
numerous and important results have been achieved on non-perturbative aspects of 
strongly interacting systems. Remarkable applications to phenomena which occur 
in statistical or condensed matter systems where the dimensionality is effectively 
reduced have been made (see for instance 0, ||). In addition, the large variety 
of two-dimensional exactly solvable models - obtained both in the scale-invariant 
regime described by Conformal Field Theories or along the integrable directions 
which depart from them - have often been the theorist's ideal playground for un- 
derstanding (at least qualitatively) phenomena which take place in high-energy 
physics. The S"-matrix approach proposed by Zamolodchikov |3| for studying the 
massive regime of integrable Renormalization Group flows has been particularly 
successful in answering important questions concerning the spectrum of QFT, the 
particle interactions and the behaviour of correlators of their local fields. 

The aim of this paper is to address the analysis of one of the most striking aspects 
of scattering experiments, i.e. the occurrence of resonances. The physical effects 
induced by resonances have been the subject of some interesting publications in the 
past. Al. Zamolodchikov, for instance, has shown in || that a single resonance state 
can induce a remarkable pattern of roaming Renormalization Group trajectories 
which have the property to pass by very closely all minimal unitary models of 
Conformal Field Theory, finally ending in the massive phase of the Ising model. 
Another model, this time with an infinite number of resonance states, has also been 
considered by A. Zamolodchikov in relation with a QFT which is characterised 
by a dynamical Z 4 symmetry ||: this is a theory with multi-channel scattering 
amplitudes which are forced by the Yang-Baxter equations to be expressed in terms 



1 



of Jacobian elliptic functions. 

The model we propose in this paper is quite close in spirit to the latter one 
by A. Zamolodchikov || but has the important advantage of being simpler and 
therefore amenable for a more detailed analysis of its off-shell properties. More 
specifically, we will consider an integrable QFT made of a stable, self-conjugated 
scalar excitation A of mass M. but with an infinite series of resonance states which 
emerge as virtual unstable particles in the scattering amplitudes. As a function 
of the rapidity variable, its two-body S-matrix is expressed in terms of Jacobian 
elliptic functions, which present a periodicity along the real axis. This causes the 
appeareance of the infinite number of resonances. In ordinary four-dimensional 
QFT, resonances are responsible for sharp peaks observed in the total cross-section 
of the scattering processes as a function of the energy or equivalently in abrupt 
changes of the phase-shifts. In two-dimensional integrable QFT, on the other 
hand, production processes are forbidden and therefore the concept itself of "total 
cross section" is not particularly usueful. We will show, however, that resonances 
are associated to rapid jumps of the phase-shift. 

The presence of an infinite number of resonance states deeply affects the physical 
properties of the system, the analytic structure of the S-matrix and in particular 
the short-distance behaviour of the correlators. This can been seen by studying the 
spectral series representations of the correlators, based on the Form Factors of the 
theory. Concerning the analysis of the Form Factors themselves, the main novelty 
consists in the infinite number of solutions of the Form Factors equations which can 
be found in this case. The infinite spectrum of operators defined by these solutions 
may be related to the irrational nature of the underlying Conformal Field Theory 
reached in the ultraviolet limit. 

The paper is organised as follows: in sect. 2 we discuss the general properties of 
our model. Its S'-matrix is the simplest example with an infinite series of resonance 
states. Since it reduces to the Sinh-Gordon S'-matrix in the limit Z — ^ (Z is the 
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module of the Jacobi elliptic functions), the QFT defined by such an S-matrix will 
be referred to as Elliptic Sinh-Gordon model (ESG). In sect. 3 we briefly outline 
the main properties of the Form Factors of our relativistic integrable QFT and 
discuss the new features which arise in a theory with double periodicity. Section 
4 is devoted to the analysis of the Form Factors of the Elliptic Ising Model, i.e. a 
model with infinite resonance states which reduces to the usual Ising model in the 
limit I — > 0. This model may be obtained by a particular analytic continuation of 
the coupling constant of the ESG model. In sect. 5 we address the computation of 
the Form Factors of the Elliptic Sinh-Gordon model: we present the computation 
of the minimal Form Factor F(f3) and investigate its analytic structure. By using 
the functional equation satisfied by F(/3), we derive the recursive equations of the 
Form Factors and obtain their first solutions. The main result of this section is that 
the solutions of the FF equations for a given number of external particles span an 
infinite dimensional vector space. In sect. 6 we present our conclusions. The paper 
also contains three appendices: appendix A presents useful mathematical identities 
used in the text, appendix B contains the derivation of the Fourier series of the S'- 
matrix and finally appendix C presents the calculation of the minimal Form Factor 
of the Elliptic Sinh-Gordon model. 

2 Elastic S'-matrix with Resonance States 

We consider a two-dimensional integrable QFT describing a stable, self-conjugated 
particle A of mass Ai. This theory is assumed to be invariant under a Z2 symmetry 
realised by A — > —A. Its on-shell properties are encoded into its elastic S'-matrix. 
In virtue of integrability, the S-matrix satisfies the factorization condition |7j], there- 
fore it is sufficient to focalise only on the two-body elastic scattering process. Let 
s — (Pi + P2) 2 be the Mandelstam variable of the scattering process and S(s) the 
two-body elastic scattering amplitude of the particle A. The function S(s) is usu- 
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ally expected to have two elastic branch-cuts along the real axis at the two-particle 
thresholds s < and s > AM 2 . These branch-cuts can be unfolded as follows. Let 
(3 be the rapidity variable which parameterises the relativistic dispersion relations 

E = Mcosh(3 ; p = Msinh(3. (2.1) 

The Mandelstam variables s and t of the two-body scattering processes are then 
expressed by 

s = ( Pl + P2 ) 2 = 2M 2 [l + coshp 12 \ ; 

t = (Pi - Vi? = 2M 2 [1 + cosh( OT - (3 12 )\ . 
(Pn = Pi — fa)- Hence the upper and lower edges of the branch-cut which starts at 
s = AM 2 are mapped onto the positive and negative real semi-axes of (3 respectively. 
The upper and lower edges of the other branch-cut are conversely mapped onto the 
positive and negative values of (3 along the line Im/5 = in. The S'-matrix S({3i 2 ) 
becomes a meromorphic function in the /?i 2 -plane. It satisfies the unitarity and 
crossing symmetry conditions 

S(P)S(-P) = 1 ; 
S(i7r -(3) = S{(3) . 

These equations automatically imply that S(f3) is a periodic function along the 
imaginary axis of the rapidity variable, i.e. 

S((3 + 2m) = S((3) . (2.4) 

This periodicity simply expresses the double-sheet structure of the Riemann sur- 
face with respect to the Mandelstam variable s where the scattering amplitude is 
meromorphic. The physical sheet can be taken then as the strip < lm(3 < tt and 
the unphysical one as the strip —in < Im/3 < 0. 

Let us now assume that S(f3) also presents a periodic behaviour along the real 
axis of the rapidity variable (3 with a period T, i.e. 

S(J3)=S(J3 + T) . (2.5) 



This equation, combined with the analyticity properties of S((3), has far-reaching 
consequences. In fact, the simultaneous validity of the two periodicity conditions 
(|2.4|) and fl2.5|) forces S({3) to be an elliptic function. Accordingly, the plane of 
the rapidity variable /3 becomes tiled in terms of the periodic cells of the function 
S(/3) (Figure l.a). On the plane of the Mandelstam variable s the S'-matrix is not 
periodic. However, the existence of a double periodicity on the /3-plane induces an 
infinite swapping between the two edges of the elastic branch-cuts on the s-plane. 
The swaps are located at the infinite set of values s n = 2Ai 2 (l + cosh(nT/2)) and 
s n = 2M 2 (l + cosh(i7r - nT/2)) for the (4A4 2 ,+oo) and (-00, 0) branch-cuts, 
respectively. In fact, take (— ^, ~) as fundamental interval of the real periodicity in 
f3 and consider, for instance, the branch-cut (4M 2 ,+oo) in the variable s. At the 
initial location s = 4M 2 ((3 = 0) the two edges join together. Starting from (3 = 
and moving onto the negative values of (3, we proceed along the lower edge of the 
branch cut, as far as (3 > — -|. When this value is crossed and (3 reaches the region 
(— T, — I"), due to the identification (3 = (3+T, we are brought back to positive values 
of (3, < (3 < ^. On the other hand, this interval for (3 parameterises the upper 
edge of the branch cut. Therefore, at (3 = — ^ a first swap of the two edges of the 
branch cut has occured. If now (3 decreases, we keep moving along the upper edge 
as far as (3 is greater than zero. When the rapidity crosses this value (now identified 
with (3 = T) a second swap occurs and the original configuration of the two edges 
is restored. The pattern described in details for < (3 < T, indeed reproduces 
itself with period T in the rapidity variable. In particular, since the point f3 = is 
identified with an infinite sequence of points (3 n = nT, the two branches return to 
the initial joined configuration an infinite number of times. For the same reason, 
the swap at (3 = —j is reproduced at every value fon+i = (2n + l)?. The situation 
is completly analogous for the branch-cut of the t-channel which originally lies on 
the real negative semi-axis of the s variable. Hence, altogether in the plane of the 
Mandelstam variable s we have the analytic structure of the S'-matrix shown in 
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Figure l.b. 

If S(/3) is an elliptic function, it must have poles and zeros, unless it is a constant. 
However, the existence of poles should not spoil the causality properties of the 
scattering theory. From a mathematical point of view, this is equivalent to requiring 
a further condition on the analytic structure of S(/3), namely that this function 
should not have poles with a real part within the physical sheet. 

It is now easy to determine the simplest iS"— matrix which satisfies all of the above 
constraints, i.e. the unitarity and crossing equations (|2.3| ), the periodicity equation 



( |2.5|) and the causality condition. It may be written as 



sn (m±^) C n (m*!2l) dn ( 
S(3a)= ± - I ± - ' ± - J (2 6) 

sn f K(/3+»ra) \ cn / KQ3-fara) \ dn ( Ktf-ina) ^ 



where sn(x), cn(x) and dn(x) are the Jacobian elliptic functions and K is the com- 
plete elliptic integral, both of modulus I (see Appendix A for a review of the main 
properties of the Jacobian elliptic functions and some of their useful identities). By 
using the addition theorems ( |A.6[ ), the above function can be equivalently expressed 
as 

sn (2Kj£) + S n(2Ka) , s 

S(/3,a)= \* R { -. 2.7 

V ' sn (») - sn(2Ka) V ; 

Exploiting some of the identities listed in Appendix A, it is easy to check the validity 
of eqs. ( |2.3| ). Moreover, the S-matrix satisfies the real periodicity condition ( |2.5| ) 
with period 

K' 

r = 7r— , (2.8) 

where K 7 is the complete elliptic integral of the the complementary modulus V = 
(1 — l 2 ) 1 / 2 . The parameter a entering the expression ( |2.6| ) may be regarded as the 
coupling constant of the model. It must be selected so as not to introduce poles in 
the S'-matrix which have a real part in the physical strip. As we will see below (eq. 
fl2.11p), this is ensured by taking a to be a real positive number in the interval [0, 1]. 



As a matter of fact, this interval can be further reduced in virtue of an additional 



property of the S-matrix ( |2.7| ). In fact, the S-matrix ( j2.7| ) is invariant under the 
replacement a — > 1 — a, 

S(P,a) = S(J3,l-a) , (2.9) 

so that a can be restricted to the interval [0, |]. For a = we have 5 = 1, i.e. a 
free theory, whereas the point a = ~ may be regarded as the self-dual point of this 
scattering theory. In the following, unless explicitly stated, a will be always taken 
as < a < 

Notice that by taking the limit I — > 0, the real period T goes to infinity and the 
S'-matrix fl2.7|) reduces to the one of the Sinh-Gordon model || 

sinh/J-isinvra 

£sh G (/M = • (2.10) 

smh p + % sin ira 

For this reason, the QFT defined by the S'-matrix ( f2.7| ) may be simply referred to 
as the Elliptic Sinh-Gordon Model (ESG). 

Let us consider the analytic structure of the S-matrix (|2.6| ), i.e. the pattern of 
its poles and zeros. Poles are either those of the numerator or come from the zeros 
of the denominator in (|2.6|) . Conversely, zeros are either those of the numerator or 
come from the poles of the denominator. Using the results listed in Table 1 and 
taking into account possible cancellations of poles against zeros, we are left with 
the following infinite set of simple poles and simple zeros 



Poles : (3 m ,n = — ina + Im-ni + nT ; 

Pm,n = ina + (2m + l)m + nT 



(2.11) 



(2.12) 



Zeros : /3 mj „ = ma + 2miii + nT ; 

An,n = — 27ra + (2m + + nT 
where m, n G Z. 

By choosing in the (3 plane the rectangle with vertices (m — ^, in + ~, —in — 
^,—in + |) as a fundamental domain, the analytic structure of the S-matrix is 
shown in Figure 2, with two zeros and two poles along the imaginary axis and 



then repeated on the other cells of the /3-plane by periodicity. Hence, the simplest 
S'-matrix with double periodicity is realised by an elliptic function of order two. 
With the above choice of the values of a, the poles of the scattering amplitude 
which have a real part in the /3-plane all lie on the unphysical sheet. Instead in the 
physical sheet there are only zeros. The poles on the unphysical sheet with a real 
part correspond to a series of resonance states JZ n with masses and decay widths 
given respectively by 



M res (n) 



n 



2M cosh (n |) cos {f) 
2M |sinh (nf) sin {jf) \ 



(2.13) 



In the limit when the modulus I of the Jacobian elliptic function goes to zero (i.e. 
T — > oo), the mass of all resonances becomes infinitely heavy and therefore these 
states completly decouple from the theory. 

It is also useful to express the S-matrix (|2.6| ) in an exponential form, a repre- 
sentation which will result particularly convenient for the computation of the Form 
Factors of the theory presented in the next sections. As shown in details in Appendix 
B, in the strip < |Im/3| < na we have the following equivalent representation of 
the S-matrix of the ESG model 



exp 



+00 1 
^— ' n 

n=l 



cosh 



nn 2 (l-2a) 
T 



COsh (2£) 



:-i) n 



sm 



2rm 



(2.14) 



If we now use the Fourier series identity (|A.13|) in order to express the factor —1 in 
front of ( gJgj) as 



— 1 = exp 



n=l v 7 



(2.15) 



the expression ( 2.14Q can be equivalently written as 



S((3) = exp 



+00 . 
^ n 



sinh 



sinh 



n(l— a)ir z 
T 



n=l 



cosh (^) 



sm 



2nir 



(2.16) 
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In the limit T 



oo, the sum in 



can be converted into an integral and 



therefore one recovers the exponential representation of the Sinh-Gordon S-matrix 



Sshc(P,a) = exp 



-4i 



dxsinh(f )sinh(f(l-a)) 



x 



cosh | 



sin 



x(3 



71 



(2.17) 



By using the expression ( |2.16|) we can read the phase-shift 5 {(3) of the model 
defined by 



sys) 



,2iS(J3) 



(2.18) 



The phase-shift is obviously a quantity which is defined modulus multiplies of ix. 
Therefore, by adding to 5(j3) a convenient multiplies of tt we can always consider 
5(P) to be a monotonic function of (3. When f3 approaches the locations of the 
resonances Re/3 n = nT, the phase-shift 8 {(3) sharply increases its value by n (see 
Figure 3). 



A particularly simple but fascinating example of integrable theory with an infi- 
nite number of resonances is given by the S'-matrix S = — 1, when expressed in the 
form ( J2.15|) . The phase-shift of such a model can be taken as the staircase function 
with values given by all multiplies of | and jumps localised in correspondence to the 
points (3 n = nT/2 (Figure 4). A way to understand this unexpectedly rich structure 
hidden in S = —1 is to notice that the scattering theory described by ( J2.15[ ) can 
be obtained as a particular analytic continuation of the elliptic expression fl2.7|) . 
Indeed, the S'-matrix ( |2.7|) approaches the value —1 every time the parameter a 
takes one of the infinite values 

T 

a c = m + (2n + l)i — , m G N, n G Z (2.19) 
2tt 

corresponding to the locations of the poles of sn(2Ka). At these critical values of 
a, the infinite number of poles and zeros of (|2.7f ) cancel each other over the entire 

9 



/3-plane, giving rise to the constant function S = — 1. In Figure 5 this pattern is 
described in details for the simplest choice a c = —iw-- In this case, for a theory 
defined at an infinitesimal displacement a = a c + 5 from the critical value a c , the 
location of poles and zeros is 

Poles : 8m li = (2n — 1)^ + Ivmxi — iixb ; 

(2.20) 

&£l = (2n + l)f + (2m + 1)tti + m5 , 
Zeros: 0mn= (2n + 1)^ + 2rrnri + mS ; 

(2-21) 

/4i = (2n - l)f + (2m + l)m - iir6 , 
where m,n G Z. As shown in Figure 5, in the limit 5-^0 cancellations occur 
between fl mn+ i, f3m}n poles and ftmji, f3 m ] l+ i zeros. The S'-matrix S = — 1 defined 
in this limit inherits all the analytic structure of the ESG model. In particular, 
the infinite number of odd values /?2n+i = (2n + l)T/2 corresponding to rapid 
increasings of the phase-shift can be interpreted as associated to resonance states 
coming from the theory ( |2.7|) defined at a = a c + 5. Near the values j3 = in + fon+i 
the ^-matrix becomes 

S(P) - — R — R , (2-22) 

P - Z7T - P 2n +1 

which corresponds to the diagram in Figure 6. In general, at the poles = iu^ 
corresponding to a bound state Ak produced by the fusion of two particles A4 and 
Aj, the S'-matrix satisfies 

-i lim W - iu%) S{P) = (r^ 2 , (2.23) 
where is the three-particle vertex on mass-shell. Therefore, comparing eqs. 



( |2.22p and ( |2.23| ), we can interpret the resonances R n as "bound states" of the 



particle A but on the second sheet, with the three-particle vertex on mass-shell 
given by = i\fn~5. The imaginary value of the three-particle vertex indicates 
that the QFT presents some non-unitary features. As bound states of the Z 2 
odd particle A, the resonance states R n are even under the Z 2 symmetry of the 
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model. Near the values (3 = /?2n+i the S-matrix behaves as in eq. ( 2.22 ) but with 
an opposite sign for the residue. These poles correspond to resonance states in the 
crossed channel. 

The even values fli n = nT corresponding to the other rapid increasings of the 
phase -shift have, on the other hand, a different physical interpretation. They do not 
come from genuine bound state resonances but simply correspond to "kinematic" 
resonances associated to the infinite replicae of the initial joint-point of the elastic 
branch-cuts on the s-plane (see Figure l.b). 



The QFT theory described by the S'-matrix ( 2.15 ) will be referred to as the El- 
liptic Ising model (EIM). The detailed study of its operator content will be discussed 
in Section 4. 



3 Form Factors in Integrable QFT 

The on-shell dynamics of a system is encoded into its S'-matrix. In order to go off- 
shell and compute the correlation functions of the model, one can use their spectral 
series representation. This means that we need to compute the matrix elements 
of the quantum fields between asymptotic states. It is sufficient to consider the 



following matrix elements, known as Form Factors (FF) [JTTJ, [12 



F n ° (&,&,..., /3 n ) = (0 I 0(0,0) | (3 u (3 2 ,...,(3 n ) in . (3.1) 

In this section we will briefly outline the main properties of the FF in a generic two- 
dimensional integrable QFT and we will comment on some features which arise in 
the case of an S-matrix with a real period. A more detailed discussion of these 
features is postponed to the next sections. For simplicity we only consider the FF 
of scalar and hermitian operators O. 

The FF of a QFT are known to be severely constrained by the relativistic sym- 
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metries. First of all, they only depend on the difference of the rapidities 

Ff(fa, /3 2 , . . . , (3 n ) = F?(/3u, Pn, .... %,•••) , i < j ■ (3.2) 



Except for pole singularities due to the one-particle intermediate states in all 
sub-channels, the form factors F n are expected to be analytic functions inside the 
physical strip < Im/?^- < ir. Moreover, they satisfy the monodromy equations 

O0 

F n (Pi, ■■, Pi, Pi+i, ■■, Pn) = F n (Pi, .., p i+ i, Pi, ..,P n )S(Pi - p i+ i) ; (3.3) 
F°(p 1 + 27ri,..,p n _ 1 ,p n ) = e 2 ™F°(p 2 ,...,p n ,p 1 ) . 

The parameter uj is the index of mutual locality of the operator O with respect to 
the particle state A (u — 0, 1/2 for local and semi-local operators, respectively). 
In the simplest case n = 2, the above equations become 

F°{0) = F°(-P)S(P), 
F°(tn-p) = e 2mw F°(iir + p) . 

By exploiting the factorization properties of integrable QFT, the general solution 
of ( |3.3|) can be found by making the ansatz 



Fn(Pi, ...,Pn)= K n (P u ...,P n )H F min (P l3 ) , (3.5) 

i<j 

where F min (P) is a solution of eqs. ( |3.4| ), analytic in < Im P < n. The remaining 
factors K n therefore satisfy the monodromy equations with 5 = 1, i.e. they are 
completely symmetric, 27ri-periodic functions of the Pi. Moreover, they contain all 
the physical poles expected in the Form Factor under investigation. 

In the case of a period S-matrix, the general solutions F n of the monodromy 
equations are not a-priori constrained to have any real periodicity in the rapidity 
variables. This can be easily seen from their parameterization ( |3.5| ); the monodromy 
equations with S = 1 which determine the functions K n do not contain any infor- 
mations about the periodic nature of the dynamics under consideration. However, 
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for the two-particle FF the solutions of eqs. ( |3.4| ) may be chosen to satisfy the extra 
condition 

F 2 °{(3 + T) = ±F 2 °{(3) . (3.6) 

In order to prove it, consider the first of eqs. ( p.4|) and make the substitution (3 — > 
(3 + T. Since S(f3) does not change under this shift in /3, if F 2 ((3) is a solution of 
eqs. ( |3.4|) , then F 2 (/3 + T) is also a solution 



F 2 °((3 + T) = S((3)F?(-(3-T) 



(3.7) 



By taking the ratio of eq. ( |3.7| ) and the first of fl3.4| ) we have 

if(/3 + T) F?{-(3-T) 



(3.1 



F°(f!) F°(-0) 

Therefore, if we assume F°(j3 + T) to be proportional to F°(/3), i.e. F°(j3 + T) = 
aF°(P), then eq. fl3.8|) implies a 2 = 1, i.e. a = ±1. In particular, for a = — 1 there 
is a doubling of the real periodicity in the two-particle Form Factor with respect 
to the S'-matrix. 

It is simple to obtain an explicit T-periodic solution of eqs. ( |3.4| ) in the case of 
a local operator (a? = 0) when the T-periodic S"-matrix is expressed in the form 

'2nn 



S({3) = exp 



— sin 

n 



n=l 



T 



(3.9) 



Denoting such a solution as F min (P), we have in fact 

1 



F min {(3) = M exp 



E 

n=l 



n 



(3.10) 



where f3 = m — (3 and M a normalization factor. 

Let us consider now the pole singularities of the FF. They give rise to a set of 
recursive equations for the F n which may be particularly important for their explicit 
determination. As functions of the rapidity differences ftij, the FF present in general 
two kinds of simple poles. The first family of singularities come from kinematical 
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poles located at = in. They correspond to the one-particle intermediate state in 
a subchannel of three-particle states which, in turn, is related to a crossing process 
of the elastic S-matrix. The corresponding residues give rise to a recursive equation 
between the n-particle and the (n + 2)-particle Form Factors [1^ 



-i lim^(/3 - 0)F° +2 + in, P, A, (3 n ) = 

(1 - e 2 ™ n"=i S (P~ Pi)) F n(Pii •••>&) ■ (3.H) 

For a theory with a real periodic S'-matrix we still expect the appearance of kine- 
matic poles located at P^ = in with residue given by the above equation, whose 
physical origin is the same as in the non-periodic case (see, for instance, the dis- 
cussion in ||1 711). However, in this case the presence of a periodic S'-matrix on the 
r.h.s. of eq. (|3.11|) gives additional constraints on the behaviour of the FF near the 
pole. Precisely, if we perform the shift Pi — > $ + T, i = 1, • • • , n in eq. ( |3.11| ), the 
r.h.s. is left invariant due to the T-periodicity of S and the traslation invariance of 
F n (the same conclusion is also obtained by shifting (3 — > (3 + T and + T on 

the l.h.s. of ( |3.11| )). Hence, in this case the residue of F® +2 in (/3 — 0) — in has to 



be a T-periodic function of the n + 1 rapidities f3, Pi, ■ ■ ■ , f3 n . Although a sufficient 
condition to implement this requirement would be F n (Pi, ...,P n ) periodic in all its 
variables, this is not a necessary condition as we will show in Sect. 5 by providing 
explicit examples. 

In F n there may also be another family of poles related to the presence of bound 
states. These poles are located at the values of Pij where two particles fuse into a 



third one. Let Pij = iu^ be one of such poles associated to the bound state Ak in 



the channel Ai x Aj and r*- the three-particle vertex on mass-shell. In this case, 
the residue equation for the S-matrix is given in ( |2.23| ) and, correspondingly, for 
F„ we obtain |T3 



-i\^eF^ +1 (P + iui k + ^P-ivf jk -^Pi,...,p n _i) = T k l3 F°{P,Pi,...,p n _i) , 

(3.12) 
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where u c ab = (n — u c ab ). As discussed in sect. 2, in the periodic case we cannot have 
genuine bound states. In fact, a bound state pole in the S'-matrix at = would 
be associated to an infinite chain of poles at = ufj + nT, n e Z with a consequent 
violation of causality. However, by considering the resonances as "bound states" 
on the second sheet, the above equation may be nevertheless used to extract some 
useful informations. 

In general, the FF of real periodic models are then expected to be ruled only by 
the recursive equations coming from the kinematical singularities, eq. (|3.11| ). In the 
particular case of the ESG model, this implies the complete decoupling between the 
two chains of FF for even or odd number of particles. This reflects the Z 2 symmetry 
of the model which can be used to classify its operators. 

Finally, we remind that in the usual non-periodic scattering theories the FF 
may be further restricted by extra conditions related to their asymptotic behaviour 



[[12] , PI and their cluster properties [10], ^ |15| . These conditions are often sufficient 
to determine uniquely the solutions of the recursive equations for each choice of 
the operator O. In the periodic it will be clear later, we deal instead with 

periodic functions for which the lack of asymptotic extra conditions leads eventually 
to an anavoidable arbitrariness in the determination of the FF. 



4 Form Factors of Ising Model Realizations 

The QFT defined by the S'-matrix S — — 1 is usually identified with the thermal 



deformation of the critical Ising model with central charge C = | fl6| . In fact, the 
solution of the FF equations in the space of hyperbolic functions permits the full 
reconstruction of all correlators of this statistical model and a precise check of its 



ultraviolet limit ]TTj, [17], |T8j . In order to point out the important differences coming 
from an elliptic interpretation of this scattering theory, it is worth summarising first 
the main results for the correlation functions of the standard thermal Ising model. 
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4.1 The Standard Thermal Ising Model 



With S — — 1, the simplest solutions of eqs. ( |3.4|) in the space of hyperbolic functions 
are given by 

f /(/?) = sinh | , if u = ; 
[ JF(/3) = tanh f , if a; = \ . 
The kinematical recursive equations for the FF are particularly simple in this theory 

- i lim - (3)F° +2 (f3 + ztt, (3, fa, ft, ... , = (l " e 2 ™ (-1)") F n (fr, ...,/?„) . 

(4.2) 

Hence, the FF of local operators (u = 0) do not have kinematical poles when n 
is an even integer whereas their residue is independent of the number of external 
particles when n is an odd integer. The situation is reversed for the FF of semi-local 
operators (u — |). 

All fields are naturally divided into two classes. The first consists of Z 2 even, 
local or Z 2 odd, semi-local fields. Local operators of this class have non-zero FF 
only between states with a finite number of external particles. The first represen- 
tative of this class is given by the trace of the stress-energy tensor 0(x). This is a 
local field, proportional to the energy operator e(x) of the thermal Ising model 

Q(x) = 2nre(x) , (4.3) 

where r = (T — T c ) is the displacement of the temperature from its critical value, 
related to the mass of the model by A4 = 2-kt. This operator has only a nonzero 
two-particle FF, given by 

F°(f3 1 -p 2 ) = 2nM 2 /(/? ;~f 2) . (4.4) 

The second class consists of the Z 2 odd, local or Z 2 even, semi-local operators, 
as for instance the magnetization operator a(x) and the disorder operator fi(x). In 
the high-temperature phase of the model, cr(x) is an odd, local operator. Its FF 
give rise to the infinite sequence F 2n+l on all odd numbers of external particles. 
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Taking into account the kinematical poles at (3^ = in, the explicit expression of 
these FF is given by 

2n+l 

= H 2n+1 ft FiPij) . (4.5) 

i<j 

H 2n +i is a normalization constant, given by H 2n +i = i n Hi, where Hi is the one- 
point FF, (0\a(0)\A) = Hl 

The disorder field n(x), in the high-temperature phase, is an even semi-local 
field. Its FF extend on an infinite sequence on all even numbers of asymptotic 
states and their expression is given by 

2n 

F 2 ~ = HtoYlFiPy) , (4.6) 

i<j 

with the normalization constant H 2n = i n H and H the vacuum expectation value 
of this operator, (0|/x(0)|0) = H . Switching between the high and low temperature 
phases, the Z 2 odd and even fields a(x) and fi(x) simply swap their role. 

The above explicit expressions of the FF allow to obtain exact informations on 
the thermal Ising model. We can compute, in particular, its correlation functions 
and related quantities. The simplest correlator is given by 

G{x) = (6(x)6(0)> , (4.7) 

(x = M.R). By inserting a complete set of intermediate states, it can be expressed 
as 

G (x) = I /'^l^| F e (/gi _ /g2) |2 e ^(cosh/ 3l+ cosh/ 32 ) = (4g) 
2 J 2tt Itx 

= M 4 [Kl{x) - Kl{x)] , 

where Ki(x) are the ordinary Bessel functions. In the ultraviolet limit x — > 0, the 
above expression reduces to 

M 2 

G(x) ~ ^ , (4.9) 

from which we can read the anomalous dimension of the energy operator in its 
conformal limit, r/ e = 1. The knowledge of G{x) also allows to determine the 
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central charge of the conformal field theory which arises in its ultraviolet limit. To 



this aim, we can use the c-theorem sum rule |20], 

C uv = - A [ d 2 x\x\ 2 (Q(x)Q(0)) 



(4.10) 



By inserting the explicit expression fl4.8|) of G(x), the above formula gives C uv = 5. 

The same set of conformal data can also be extracted by looking at the ground- 
state energy of the theory on an infinite cylinder of width R. The exact expression 
is obtained by means of the Thermodynamic Bethe Ansatz (TBA) |22| and for the 
effective central charge 

C = C 24A rri j T1 , 



we haveQ 



C(x) 



1 

2 




, 1 1 i 

m - + - + In 7r - 7s 

x 2 



(4.11) 



(4.12) 



>-2n-l 



n + 1 



)£(2n + l) - 



X' 



2 \ n 



7T Z 



The power series of the above expression should match with its evaluation in terms 
of Conformal Perturbation Theory given by 



Cper (■£) 



(4.13) 



71=1 



where 



12- 



-l) n 



nl 



R 2(l-n)+n Ve / ^ A 



n-1 



cyl 



:(0) JJe(^)d 2 Wi |A min )) c 

3=1 



(4.14) 



and ((•••)) denotes the conformal correlators on the cylinder. 

Since for the critical Ising model A min = and the energy operator e(r) is odd 

under duality, we have c n = for odd n. Hence, eq. (|4.12 ) fixes the central charge 

to be C uv = ^. Moreover, from the comparison of the power series in (gjj) with the 
*In this equation je is the Euler-Mascheroni constant and £(s) the Riemann zeta function. 
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perturbative expansion ( 4.13 ), we have i] t = 1, in agreement with the determination 
of this quantity obtained by the short-distance limit of the correlator ( |4.8|) . 

Let us now consider the correlation functions of the Z 2 odd, local and even, 
semi-local fields, 

G(x) = (a(R)a(0)) ; 

G+(x) = (M^)Mo)> , 

given by 

G ± M = E^/f -%F^e-^f», (4.16) 

n 

where the sum runs on odd numbers for G + (x) and on even ones for G_(x). As 
shown in [|17], 19| , eq. ( |4. 16|) provides an integral representation of a particular so- 
lution of the Painleve' equation. An exact analysis of this solution allows, in par- 
ticular, to study its short-distance limit x — > with the result 

„ , . const , A 
G±{x) ~ — - . (4.17) 
R* 

The above expression determines the anomalous dimensions of the Z 2 odd fields to 
be T\ a = T) lk = §. 

The anomalous dimension of these fields can be equivalently extracted by means 
of the 7] sum-rule fT4] 

r]rS> = -—L^ f d 2 x (6(x)$(0)) . (4.18) 



2vr($) 

This formula can be applied in the high-temperature phase to obtain the anomalous 
dimension of the disorder field and in the low-temperature to obtain the anomalous 
dimension of the magnetization operator. By inserting a complete set of states 
between the two operators and noticing that the FF of 0(x) is different from zero 
only on the two-particle state, we have 
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It is worth mentioning that the anomalous dimension of these fields can be 
also obtained by exploiting the formal analogy of expression ( |4.16|) with a grand- 
canonical partition function S(£, L) of a fictitious one-dimensional gas in a box of 



length L ~ log and fugacity z = ^ 



1/, 18]. Hence, by the state equation of 



this gas 



Y \p( z ) 

alr] 



(4.20) 



the anomalous dimension is nothing but its pressure p(z) at z = 5-. In the nearest 
neighborhood approximation of "E{z, L), the pressure p can be obtained as a solution 
of the integral equation 



2tt 



dze 



-p z 



(4.21) 



In the thermal Ising model, the above formula provides an excellent approximation 



of 77, i.e. r] a ~ 0.12529 [18 



4.2 Form Factors for the Elliptic Ising Model 

Let us discuss now the Form Factors of the theory described by the 5-matrix 
S = — 1 but this time regarded as a periodic function of (3. As it is clear from 
eq. ( |A.13| ), there is a complete arbitrariness in the choice of the periodic realization 



of — 1. Instead of working with the S'-matrix ( 2.15 ) of period T, it is more convenient 
to consider the following 2T periodic realization 



— 1 = exp 



n=l 



2n 



sin 



n(2n - l)f3 
T 



(4.22) 



which leads to a simplification in the expressions of the Form Factors of the model, 
without altering its physical content []. 

t In fact, one can always convert the expressions of the FF relative to the period 2T to those 
relative to period T by multiplying them for even, 2ni periodic functions which enter the general 



solution of (4.23). For instance, the function sn( i^jr/3) , solution of (4.23) with period 2T, can be 
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In what follows we will not address the problem of finding the most general 
solutions of all FF equations relative to this theory. We will rather concentrate 
our attention on particular classes of solutions which have the minimal analytic 
structure compatible with all the constraints. 

Let us consider first the minimal two-particle Form Factors, solutions of the 
equations 



They are determined up to arbitrary even, 27rz-periodic and 2T (anti)-periodic, 
analytic functions of j3. This time we are looking for the simplest solution of these 
equations in the space of elliptic functions. Using the properties of the Jacobian 
elliptic functions listed in Appendix A, it is easy to see that they are given by 



These expressions represent the simplest generalization of the standard results ( |4.1|) 
with no poles on the real axis and which accounts for the real periodicity of the 
system under consideration, with the property F min (P + 2T) = — -F min (/3). Their 
analytic structure consists of the following set of simple zeros and poles: 



n-P) ; 



F(m - (3) = e 2 ™ F(m + 0) . 



(4.23) 




(4.24) 



transformed by means of the even, 27ri periodic function I sn' 




+ 1 into the T-periodic 



solution 




of the same equation. 



Zsn 2 (i£ji/3)+l 
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HP) 



Zeros : f3 m ,n = 2mni + 2nT 
Poles : /3 m>n = (2m + l)m + (2n + l)T ; 

(4.25) 

Zeros : f3 m ,n = 2miri + 2nT ; 
Poles: fl m ,n — {2m + Y)Tii + 2nT , 

with m, n e Z. Notice that /(/3) has an infinite sequence of poles at (3 = in + (2n + 
1)T, located just at the edge of the second sheet. They correspond to the infinite 
sequence of the one-particle intermediate states given by the resonances R n (see 
the discussion at the end of Section 2). On the contrary, the poles of located 
at f3 = m + 2nT do not signal the presence of resonances but are simply related to 
the infinitely repeated pattern of the branch cut of the t-channel of the S'-matrix. 

Concerning the recursive equations satisfied by the FF, they are given by the 
same equation ( (4.2|) and therefore, as before, the operators fall into two different 
sectors. 

Let us consider first the Form Factors of the local Z 2 even or semi-local odd 
operators, i.e. the ones which do not have poles at fiij = in. Their general expression 
can be written in terms of the function /(/?) and they are nonzero only on a finite 
number of external states. Let us consider, in particular, the FF of the trace 
operator. Even in the periodic case we assume this operator to be related to the 
energy field by eq. ( |4.3| ). Its only non-vanishing FF is the one on the two-particle 
state which we assume to be 

F e (A-/? 2 ) = <0|e(0)|Aft> = 2tt7W 2 /( ^ ~f 2) . (4.26) 

For the two-point function of this operator we have 



G(x) = (9(i2)e(0)> = (4.27) 

_ 1 r d£id^2 \p0(g _ g \| 2 ^-x(cosh/?i+cosh/3a) 
2 J 2vr 2vr 1 K ^ ^' 
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By using the variables y± = \{@x ± P2), it can be written as 

4M 4 f +co 

G(X) = [7MF Jo ^-l^( 2 ^-)| 2 ^( 2xcosh ^) ' ( 4 - 28 ) 

Plots of this correlation function for different values of the modulus / are shown in 
Figure 7. 

Let us discuss in more detail the large and the short distances behaviour of this 
correlator. As evident from Figure 7, the large distance behaviour of the correlator 



( f4.28| ) is essentially the same for all value of the period T and is ruled by the 
mass gap A4 of the stable particle. This is confirmed by the analytic estimate of 
G(x) for large values of x: the only parameter which contains informations on the 
periodic structure of the theory (in terms of parameters of the elliptic functions) is 
the normalization in front of the exponential decay 

2 o p~ 2x 

G{x)~-(M 2 l'K) -^ + 0(e~ ix ) . (4.29) 

7T OC 

The situation is drastically different in the ultraviolet regime. While in the non- 
periodic case (I = 0) the two-particle FF is expressed by the unbounded function 
sinh |, for I 7^ it is instead a periodic, regular and limited function on the real 
axis which can always be bounded by a constant B, leading to0 

p+00 

G(x) < AM A B 2 / dyK (2xcoshy) = 2 (M 2 B) 2 K 2 {x) . (4.30) 
Jo 

This inequality implies that for / 7^ G(x) cannot have a power law singularity 
in its ultraviolet limit as, instead, it happens in the standard Ising Model (see eq. 
(fOD ) • Hence the conclusion is simply that the anomalous dimension of the energy 
operator changes discontinously from the value r\ e — 1 (for / = 0) to r] e = for all 
other finite values of the module! 

This result seems to give rise to an apparent paradox. In fact, one could argue 
that the anomalous dimension of the energy operator can be equivalently extracted 
* To obtain this result we have used eq. (6. 663.1) of p3|. 
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in terms of the TBA equations by looking at the perturbative part of the ground 
state energy. Since this quantity is still given by the same eq. (|4.12j) , it appears 
that from the TBA approach one would rather reach the previous conclusion, i.e. 
rj e — 1. How can we conciliate this apparent mismatch of the values of rj t obtained 
by the two different methods? 

The solution to this puzzle is the following. First of all, notice that the TBA 
does not compute the central charge of the theory but rather provides the evaluation 
of the effective central charge of the theory, C = C — 24A min . To compute directly 
the central charge, we have instead to employ the c-theorem sum-rule ( |4.10|) . The 
result for the central charge of the Elliptic Ising Model is reported in Figure 8. This 
figure shows that, by varying the period, the central charge of the model is no longer 
C — \ but, on the contrary, takes any value between and |, depending on the 
module I. On the other hand, since the effective central charge C is always equal 
to | for any value of the module I, there should be in the periodic realization of the 
Ising model an operator with conformal weight A m ; n < 0. Hence, the underlying 
conformal theory of the Elliptic Ising Model is no longer unitary (and rational) 
and consequently we have no reasons to argue that the perturbative coefficients c n 
vanish for n odd, as in the Ising Model. Therefore, the matching of the Conformal 
Perturbation Theory in the presence of a field with A m i n ^ with the corresponding 
perturbative part of (|4.12 ) gives rj e = 0, in agreement with the determination done 
by looking at the ultraviolet limit of the correlator. For the limiting value 1 = 1, 
we have A min = -i. 

As already discussed, the FF of even, local operators inherit from the S'-matrix 
the pole structure associated to the resonance states. These states are even under 
the Z 2 symmetry of the model and therefore they couple directly to the operator 0. 
We can extract the matrix elements (0\Q(0)\R n ) by formally applying the bound 



state residue equation (|3.12|) , with T A \ = iynS (8 — > 0). As a result we have 

(-l) n -^ = ^(010(0)1^) . (4.31) 
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Hence, (0|O(0)|i? n ) diverges as l/y/6 but its product with the three-particle vertex 
of the resonances is constant in the limit 5^0. 

Let us consider now the FF of the Z2 odd, local and even, semi-local fields, 
concentrating our attention on the magnetization and disorder operators. The 
magnetization a(x) is local and has FF on all odd numbers of external particles. A 
close solution of the recursive equations ( }4.2|) is given by 

2±\ {lT (n + i) Hi Y[ m .) , (4.32) 

where #1 = (0\a(0)\A). 

The disorder operator n(x) is semi-local and has FF on all even numbers of 
external particles 

-^J (I'fHoUmi) , (4-33) 
' i<j 

where H = (0|/i(0)|0). All the FF ( ggg ) and flOS) ) possess the kinematic pole at 
Pij = in as well as at their periodic repeated positions = ire + 2nT. 



For the two-point correlation functions of these fields we have 



with 



_(*) = 


(a(R)a(0)) 


+( x ) = 




fdf3 1 


d(3 n | ± 


J 2tt 


2vr 1 n 



(4.34) 



^)=£^ I ^■■■^|i£|V*^* (435) 



where the sum runs on odd numbers for £+(x) and on even ones for £_(x). 

As for the energy field, at large distances these correlation functions have an 
exponential decay, ruled by the mass gap M. of the stable part of the spectrum 



g-(x) - H "y^ e ' x + ( e ~ 3x "> ' ( 4 - 36 ) 

g + (x)^^(l'^H ) 2e -^ + O(e-^) . 
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The only informations on the periodic structure of the theory may enter the normal- 
ization in front of the exponentials. Let us now analyse the short-distance behaviour 
of G±(x). In absence of an exact resummation of the series ( |4.35| ), we have to rely 
on other methods for computing the anomalous dimensions of the corresponding 
fields. It is easy to see that, in contradistinction to the anomalous dimension rj e of 
the energy operator which jumps discontinuosly from 1 to by introducing a finite 
period T, the anomalous dimension r] a (or equivalently i]^) depends continuously 
on T. In fact, in the ultraviolet limit we can apply the thermodynamics approach 
as in eq. ( |4.20|) and write the following representation for the correlator 

/ x y(i) 

S = S + + g -^{MR) • (437 > 
where the pressure of the one-dimensional fictitious gas is given in eq. ( |4.21| ) . Since 
in this case |JF(/?)| 2 is a regular and bounded function along the real axis for any 
value of the module I, it follows that p{l) is a continous function of the module. 
Another way to reach the same conclusions is to analyse the exact expression for 
the anomalous dimensions as obtained by the i] sum-rule 

*<0=<fc<D = S y c cosh ,p Ha J^-y (4.38) 

The relative graph of this function is reported in Figure 9. The curve starts from 
r) — g at I — and goes then continously to zero at I — 1. 



5 Form Factors for the Elliptic Sinh— Gordon Model 

Let us now consider the quantum field theory described by the elliptic S"-matrix 
( |2.6| ) and address in its full generality the problem of computing exactly its FF. We 
will see that the FF equations admit in this case a large class of solutions which for 
I —>■ all collapse to the same solutions of the Sinh-Gordon model. As a matter of 
fact, due to the absence of constraints on the asymptotic behaviour of FF, the final 
set of solutions has indeed infinite dimension. 
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We concentrate only on the case of FF of local operators. As usual, the first step 
consists in the evaluation of the two-particle Form Factor. For any local operator, 
this is a solution of the monodromy eqs. ( |3.4|) with uj = 0. Using the exponential 
representation ( |2.16| ) for the S-matrix and the formula ( |3.10| ) we have 



F min ((3) = Afexp 



oo . 

^ n 



sinh 



sinh 



T / - 2 / 7T 



71=1 



, , sin \n—8 
cosh (of) sinh (2^) V T 



(5.1) 



where f3 = iir — (3 and M is a normalization constant. This solution satisfies 
the periodicity condition F min (/3 + T) = F min ((3). The previous expression for 
-^min can be further elaborated in order to get a final form, more suitable for the 
determination of its analytic structure and for explicit calculations. As described in 
details in Appendix C, it can be rewritten, up to a constant, as an infinite product 
of trigonometric functions 



(5.2) 



k=2 



where 

(l - 2q 6 cos (2f p) + qt) (l - 2q 5 cos (2f/?) + qfj (l - 2g 4 cos (2*/?) + q\ 
(\ - 2q x cos (2%p\ + q{A (l - 2q 2 cos fo^f) + q^j (l - 2q 3 cos (2%$} + gf 

and qi = qi(k), i = 1, • • • , 6 are the factors listed in eq. ( |C.7| ). Using the identity 
( [A. 18 ) we can also express the previous result as an infinite product of gamma 
functions 

+00 +00 

F mi M= n n y kM , (5.3) 



where 



Y k , n (l3) 



n=—oo k=0 



2 ~ t ~ 2tt 



r(* + § + ^)r(* + f-f + £)r(* + i + § + £ 



and (3 n = (3 + nT. 
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(5.4) 



The solution (|5.2[ ) is analytic in the /3-plane except for an infinite set of poles 
(see Appendix C for details) 

Poles: /3 n (k) = — ma + (2k — l)m + nT ; 

(3 n (k) = -ma + (2-k)2m + nT ; 

Pn{k) = ma + (3 — 2k)m + nT ; 

/^n(^) = i^a + [k — l)2m + nT ; 

Pn{k) = 2km + nT ; 

/5 n (fc) = (1 - k)2m + nT , 

and it vanishes at the following locations 

Zeros : (3 n {k) = — ma + 2km + nT ; 

Pn{k) = — ma + (3 — 2k)in + nT ; 

Pn{k) = ma + (1 — k)2m + nT ; 

Pn{k) = ma + (2fc — l)m + nT ; 

(3 n (k) — (k- l)2m + nT ; 

(3 n (k) = (2-k)2m + nT . 

For given /c, poles and zeros have multiplicity (k — 1) and they are drawn in Figure 
10. The minimal solution T min does not have singularities in the physical sheet but 
only single zeros in (3 = nT. The poles at (3 = —ma + nT on the unphysical sheet of 
T m m(/5) may be regarded as "bound state" poles of the FF due to the resonances. 
The solution (15.21) satisfies the functional relation 



(5.5) 



sn (2i—8) 

^ + ») F mM = 311(2if;3) i _ T 3 : (2Ka) • (^) 

This can be easily proven by studying the distribution of poles and zeros for the 
product on the l.h.s., as shown in Figure 11. When -F min (/3) is multiplied by the 
same function evaluated at the shifted point (f3 + in) a partial cancellation of poles 
against zeros occurs, leaving only a distribution of simple poles and zeros (third line 
of Figure 11). This distribution indeed corresponds to the pole singularities and 

28 



zeros of the elliptic function on the r.h.s. of eq. ( |5.6| ). Notice that we have used the 
functional equation ( |5.6|) as it stands in order to fix a convenient normalization TV 
of F min ((3). 

In order to simplify the following equations, it is convenient to introduce the 
notation 



u = i—/3 , Uij = i—pij ■ (5.7) 



rp P 1 U «j — 1 rp 

The general structure of rz-particle Form Factors is given in eq. fl3.5p in terms 
of the minimal solution ( |5.2j ). Since the theory has no bound states, the only 
poles present are the kinematic poles associated to every three-particle subchannels 
(Pij = wr). This suggests to parameterise the functions K n entering the general 
expression of the FF (|3.5|) as follows 



K n ((3 ir --,f3 n ) = H n Q n {pi, •••,/?„) TT —— , (5-8) 

- LJ - cnuj,' 

where H n are normalization constants. Notice that in the above expression the 
presence of the kinematic pole has been made explicit by the introduction of the 
periodic elliptic functions cnu^. The choice of this parameterization is suggested 
by the general requirement that the residue of F n at the kinematic pole has to be a 
periodic function of the other (n — 2) rapidity variables (see the general discussion 
in Sect. 3). Therefore, the kinematic pole at /3y = in appears toghether with an 
infinite sequence of poles located at = m + 2nT. Even though the function cnu 
is 2T periodic, the actual T-periodicity of the residue of F n (fii, ...f3 n ) at = in 
will be restored by the consequent expressions of Q n (Pi, ■■■Pn) obtained by using 
this parameterization. 

The function cn[u(/3)] satisfies the identity 

cn[u(P + ivr)] cn[u(/3)l = I 1 sn ^ cnu . (5.9) 

dnu 

The terms Q n are symmetric, 2ni periodic functions of the /3j's, which have to 
be determined by the recursive equations ( |3.11 ) satisfied by the FF. By using the 
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expression (277) for the S'-matrix and exploiting the functional relations ( |5.6| ) and 



(|5.9|) together with eq. ( |A.10|) , we obtain the following recursive equations for the 
functions Q n 

Q n +2(P + in,P,Pi,-'-,Pn) = V n (P\[3 u ---,p n )Q n ([3 u ---,p n ) (5.10) 

where 

v n {p\P,-M = ^^y* (5.11) 

Iir =1 [sn[2(u - uj] - sn (2Ka)] - QlLi [sn[2(u - uQ] + sn (2Ka)] 
^[l + dnpCu-Ui)]] 
Concerning the normalization constants if n , they have been conveniently chosen as 

V 2 / V vrF min (m) / 

V2y V vrF min (z7r) / 
with the constants H , Hi determined by the vacuum and the one-particle FF of 
the operator O. 

Let us discuss the general properties of the recursive equations ( |5.10| ). It is easy 



to see that the solutions Q n with n > 1 span in this case an infinite dimensional 
vector space. In fact, if Q n (Pi, ■ ■■,P n ) is a non-trivial solution of the recursive eq. 
( |5.10|) , the general expression 



Q n (Px,...,Pn) = Qn(Pl,-,Pn) + MPl, ■■; Pn) Wfa, Pn) (5.14) 

is still a solution if W(Pi, ...,/3 n ) belongs to the kernel of the recursive equation 

W n {P + iir,P,p 3 ,...,P n ) = , (5.15) 

i.e. W n (Pi,..., P n ) = Yl i<:j cn Ujj, and A n (Pi, ■ ■ ■ , p n ) is a generic 2iri periodic, sym- 
metric function with possible pole singularities outside the physical strip and which 
reduces to a constant in the limit I — > - the last requirement in order to reproduce 
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the standard result for the Sinh-Gordon model|. Due to the periodic nature of the 
theory under consideration, no extra constraints can be imposed on the function 
A n which then spans an infinite space of solutions. Moreover, as we will see from 
explicit examples below, there exist in general several different non-trivial solutions 
Q n which however reduce to the same solution of the FF equations of the Sinh- 
Gordon model in the limit I — > 0. In other terms, given a solution Q n of the FF 
equations of the Sinh-Gordon model, there are in general many different ways to 
extend this function for I 7^ and obtain then distinct solutions of the FF equation 
of the Elliptic Sinh-Gordon model. 

Let us concentrate now on the construction of particular solutions Q n of the 
recursive equations. 

For n — 1, Qi is a constant, as required by relativistic invariance. By absorbing 
this constant into the normalization Hi, we can take Q\ = 1. 
For n = 2a solution is given by 

QziPufo) = cnu 12 . (5.16) 

In fact, the two-particle FF of any local operator cannot have a pole at (3^ = in 

so that Qi{Pi,j32) must precisely cancel the pole introduced by the function cnui2 

in the parameterization ( ]5.8|) . For / — > this solution reduces to the corresponding 

one for the Sinh-Gordon model || [K]]. However, this is not the only acceptable 

solution with this property. We can for instance divide the function (|5.16|) by dn U12 

and obtain another admissible two-particle FF, with no pole at f3\2 = in and which 

reduces to the same function of the Sinh-Gordon model in the limit I —>■ 0. This new 

solution differs from the previous one by the distribution of zeros and poles on the 

complex /3-plane. In particular, it possesses poles at j3u = (2m + l)in + {2n + 1)T 

and zeros at (3\2 = limn + (2n + 1)T, which are pushed to infinity for / — > 0. 

Hn the Sinh-Gordon case this a-priori infinite set of solutions is severely restricted by the 
constraints on the asymptotic behaviour of F n which eventually force the function A n to be a 
constant |Tc|| . 
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Along these lines, it is easy to see that the most general solution to the recursive 
equations for n = 2 falls into the class ( |5.14| ), i.e. it may be written as 



QziPufo) = cnuis q(/3i2,l) 



(5.17) 



where q{f3u, I) is any 2ni periodic, even function with no pole at /3i 2 = iir and such 
that lim/^o — constant. Different solutions will have different distributions 
of zeros and poles but, in the absence of extra requirements on the FF of this theory, 
they are all acceptable, as far as they do not possess dangerous poles on the physical 
strip. 

Let us now discuss the three-particle FF. The function Q 3 (/?i, /3 2 , /%) must be a 
solution of the equations (|5.10|) specialized to the case n = 3 

2 



Q 3 (P + m,P,(h) 



(5.1* 



1 + dn[2(u - u 3 )] 

where we have used Qi — 1, absorbing the constant in Hi. It is easy to find two 
different non-trivial solutions of this equation with no poles in the physical sheet 
of each two-particle subcluster. The first one is given by 



+ 



cn(u i2 + u 13 ) cnu 23 cn(u 23 + u 2i ) cnui 3 
dn(ui2 + Ui 3 ) dn u 23 dn(u 23 + u 2 i) dn u i3 
cn(ui 3 + u 23 ) cnui2 



dn(ui 3 + u 23 ) dnui 2 _ 



(5.19) 



The second solution is given by 



cn 2 U12 cn 2 ui 3 cn 2 u 23 



(5.20) 



_ dn 2 U12 dn 2 u i3 dn 2 u 23 
Both expressions for T — > oo fall into the general Q 3 solution for the Sinh-Gordon 
model which is given by |T0|| . 



Ql h (Pi, fh, fa) = AiicoahPu + cosh/5i 3 + cosh/3 23 ) + M 



(5.21) 



with A 2 — Ai = 1. We notice that both the solutions are periodic in each variable, 
with period 2T, even if they satisfy the condition for the residue at the kinematic 
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pole to be T-periodic, as it follows from the recursive equation ( |3.11 ). However, in 
general nothing prevents one from finding also T-periodic solutions. 

As in the n = 2 case we can generate a pletora of solutions by multiplying the 
solutions (|5.19| ) and ( |5.20| ) by any function q(fa, fa, fa, I) which is 2ni periodic, 
symmetric, with q(/3 + in, fa fa, I) = constant and such that it reduces to a constant 
in the / — > limit. Such a function may be written, for instance, in terms of 
arbitrary powers of the expression qo(fa, fa, fa) = dnui2 dnui 3 dnu2 3 . 

As a final example, we provide explicit solutions of the recursive equations (|5.10Q 
for n = 4. In this case, we need to find Q 4 satisfying 

sn[2(u — u 3 )] + sn[2(u — u 4 )] 



Qi(f3 + m,fafa,fa) 



-Ai 



[1 + dn[2(u - u 3 )]] [1 + dn[2(u - u 4 )]] 



x Q2(fa,fa 



(5.22) 



Using some of the identities listed in Appendix A and inserting the explicit expres- 
sion for Q2 as given in eq. ( 5.16Q , the previous equation can be rewritten as 

sn(2u — u 3 — u 4 ) cn 2 (u 3 — u 4 ) dn(u 3 — u 4 ) 



—at 



Q 4 ((3 + in,f3,fa,fa) 
If we now define the functions 



[dn(2u — u 3 — u 4 ) + dn(u 3 — u 4 )]" 



Di = J^dnuy , 

i<j 



(5.23) 



(5.24) 



(5.25) 



dn(ui2 + u 34 ) dn(ui2 - u 34 ) dn(u i4 + u 23 ) , 

it is easy to determine various solutions with no poles on the physical strip. Here 
we list two of them. 
The first solution is 



Q?(fa,fa,fa,fa) 



D\D 2 



At 



V D\ D 2 + (I') 



x 



(cn U12 • cn u 34 + cn Ui 3 ■ cn U2 4 + cn u 44 ■ cn U2 3 ) x 



(5.26) 



1 cn(u i2 + u 34 ) cn(ui2 - u 34 ) cn(u i3 + u 23 ) cn(u i3 - u 23 ) 



+ 



dn(ui2 + u 34 ) dn(ui2 - u 34 ) dn(u i3 + u 23 ) dn(ui 3 - u 23 ) 
cn(u i4 + u 24 ) cn(u i4 - u 24 ) 



dn(u i4 + u 24 ) dn(u i4 - u 24 ) 
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In the limit T — > oo it reduces to the standard Sinh-Gordon solution 

ysv.., > , , , , • / . i_ P12 flu , fll3 fl\ 



[2 + cosh /3i2 + cosh fli 3 + cosh flu + cosh /5 23 + cosh /5 24 + cosh /9 3 4 



■ I u u ^34 . , Pl3 , P24 . , PU , 

z cosh — cosh h cosh — cosh h cosh — cosh 

V 2 2 2 2 2 



The Sinh-Gordon FF relative to this solution goes asymptotically to a constant in 
each variable 

The second solution reads 

D\D\ 



Qi 2) (/3i> = 
"cn(ui2 + u 34 ) 



A3 



X 



( J D? J D 2 + (/') 6 ) 2 
cn(ui 3 - u 24 ) cn(u M + u 23 ) 



dn(ui2 + u 34 ) dn(ui2 - U24) dn(ui 4 + U23) 
cn(ui2 + u 34 ) cn(ui 3 - u 24 ) cn(u M + u 23 ) 



+ 



(5.28) 



dn(ui2 + u 34 ) dn(ui2 - u 24 ) dn(u M + u 23 ) 
and in the Sinh-Gordon limit it reduces to the second standard solution of this 
model [Kj 



Q 4 5/l (A, &,&,&) 



-2, I cosh ■ cosh . cosh W» + M 



+ cosh^±M + cosh-^ + ^ 



+ cosh 



(flu + /3 23 ; 



(5. 



2 2 2 

The Sinh-Gordon FF relative to this solution goes asymptotically to zero in each 
variable fli. 

Starting from a different expression for Q 2 , one can easily generate other solu- 
tions for the four-particle Form Factors which however reduce either to ( |5.27| ) or 
( gTgp in the limit 1^0. 

Concerning the computation of higher FF, it can be performed along the same 
lines of the previous examples. As shown by the first few cases, the periodic nature 
of the theory will manifest itself in the appearance of infinite sets of solutions 
which all collapse to the standard solutions for the Sinh-Gordon model in the 
I — ► limit. As far as the modulus I is kept away from zero, different solutions 
for the n-particle Form Factors are characterized by a different analytic structure. 
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It would be interesting to understand how to discriminate among different FF on 
the basis of a detailed investigation of the zeros and poles distributions of matrix 
elements of the various operators present in the theory. Moreover, it would be 
highly desirable to have a concise and closed formula at least for a sequence of Form 



Factors F n (f3\ } ...,/3 n ) as in the Sinh-Gordon model [10|, but the determination of 
this formula has presently eluded our attempts. 

Finally, let us briefly discuss some aspects of the ESG model in its ultraviolet 
regime. Among the solutions of the FF equation we expect to find those ones which 
identify with the matrix elements of the trace Q(x) of the stress-energy tensor. This 
is a Z 2 even field, with non-vanishing FF for all 2n numbers of external particles. 
The FF of the field Q(x) can be used to compute the central charge of the model 
in its ultraviolet limit by means of the c-theorem sum rule 



oo 



C(T) = / g^O) , (5.30) 
Jo 

where Ci(/x) is given by 

x 5(y~^ j m sinh m cosh (3j — /i) . 

i i 

Assuming the uniform convergence of this series (so that we can interchange the 
sum on the index n with the integral in ( |5.30| )), it is immediate to conclude that 
the central charge of the model can take any value in the interval (0, 1) depending 
on the period T of the S— matrix. In fact, for T — > oo we recover the sum rule 
series of the Sinh-Gordon model, i.e. C(oo) = 1. Viceversa for T — > each term 
of the series goes to zero and therefore (7(0) = 0. The series (|5.31|) is a continous 



and monotonic decreasing function of T so that by varying this parameter we can 
reach any value in the interval (0,1). Therefore, in the ultraviolet limit we have 
generically a non-unitary irrational Conformal Field Theory. As a consequence, it 
seems reasonable that the existence of an infinite dimensional space of solutions for 
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the FF of the ESG theory can be traced back to the non-unitary irrational nature of 
the CFT in the UV region and to its relative infinite number of primary operators. 

6 Conclusions 

In this paper we have analysed certain massive integrable Z 2 invariant quantum 
field theories with an infinite tower of resonance states. These states are associated 
to an unlimited sequence of poles on the unphysical sheet of the S-matrix for a 
fundamental particle A and may be regarded as unstable bound states thereof. As 
a function of the energy, the phase-shift of the scattering amplitude presents sharp 
jumps in correspondence to the masses of the resonances. The S'-matrix of these 
theories has a real periodic behaviour in the variable (5. 

The Form Factors of such theories may be computed in principle along the same 
lines of other integrable relativistic models, i.e. by solving the monodromy and 
recursive equations coming from the residue condition on the kinematical poles. 
However, new features have emerged from the periodic nature of the S'-matrix: on 
the one hand, there is a very severe constraint on the residues of the Form Factors 
F n (Pi, ...p n ) at their kinematical poles = in, which have to be periodic expres- 
sions of the remaining rapidity variables (5^. Conversely, an infinite proliferation of 
solutions of the FF equations occurs, due to the lack of very stringent conditions 
on their analytic structure or to the impossibility of enforcing a given asymptotic 
behaviour in each rapidity variable 

While it is an interesting open problem to develop further theoretical criteria 
to identify the operators associated to these solutions, nevertheless the presence of 
an infinite number of them seems compatible with an other aspect of these models 
related to their behavior in the ultraviolet regime. In fact, the Conformal Field 
Theories which are approached in the ultraviolet limit are generally non-unitary 
and irrational, therefore with an infinite number of primary fields. The value of 
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the central charge may be fixed by varying the spectrum of the resonances, i.e. by 
changing the period T of the S-matrix. 

Since these QFT with resonance states may be thought of as theories with 
additional mass scales compared to those given by the spectrum of stable particles, 
it is an interesting general problem to determine how these extra mass scales affect 
the different regimes of the theory. From our analysis it has emerged that the large 
distance behaviour of the model is hardly affected by them, being essentially ruled 
by the mass gaps of the stable part of the spectrum. On the contrary, the ultraviolet 
properties of the theory are deeply influenced by the infinite tower of unstable states 
which, in the ultraviolet regime, tend to smooth the short-distance divergencies 
and to decrease their related quantities. As a consequence, ultraviolet singularities 
are in general less severe, and the central charge of the underlying Conformal Field 
Theory takes lower (in general irrational) values than in the absence of the resonance 
states. It would be interesting to further investigate the ultraviolet properties of 
these theories with resonances by means of the Thermodynamics Bethe Ansatz. 

More generally, it would be highly desirable to understand whether these theories 
allow for a Lagrangian description, and to find possible applications to statistical 
mechanics systems. Moreover, the possibility to define other periodic models cer- 
tainly deserves more investigation as well as the possible computation of FF for the 
Z 4 -model of Ref. f§. 
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Appendix A 



In this appendix we collect some useful formulas relative to the Jacobian elliptic 
functions as well as mathematical identities used in the text. The reader may 



consult |24| for further details. 

The Jacobian elliptic functions are double-periodic functions and have two sim- 
ple poles and two simple zeros in a period parallelogram. Let 

k(o=/ ? . d :. 2 , (a.i) 

Jo v 1 — I 2 sin a 

be the complete elliptic integral of modulus I and K'(/) = K(Z') the complete 
elliptic integral of the complementary modulus with I 2 + V 2 = 1. Setting q = 
exp [— 7rK'/K], the definition of the Jacobian elliptic functions snu, cnu and dnu is 
given by 

2q3 . -KU 1 - 2q 2n cos ^ + q 4n 

SI1U ~ 7f Sm 2K l = ll_2q2n-i C os^ + q 4n - 2 ' 

2Wq^ ttu ^ l + 2q 2n cosf + q 4n 
CI1U = ~7T C ° S 2K II l-2q^cos^ + q^ ' ^ 



dnu 



n=l K 

1 + 2q 2n - 1 cos ^ + q 4n " 2 
2q 2n - 1 cos ^ + q 4n " 2 



n= i " K 
The parity properties of the Jacobian elliptic functions are 

cn(— u) = cnu ; sn(— u) = — snu ; dn(— u) = dnu , (A. 3) 

whereas their value in u = is 

cn(0) = dn(0) = 1 ; sn(0) = . (A.4) 
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The period, zeros and poles of these functions are summarised in the following table 



Function 


Periods 


Zeros 


Poles 


snu 


4mK + 2nK'i 


2mK + 2nK'i 


2mK + (2ra + l)K'i 


cnu 


4mK + 2n(K + K'i) 


(2m + 1)K + 2nK'i 


2mK + (2n + l)K'i 


dnu 


2mK + AnK'i 


(2m + 1)K + (2n + l)K'i 


2mK + (2n + l)K'i 



Table 1 

The change of argument of these functions are ruled by the table 



u* = u + K 


u + iK' 


u + 2K 


u + 2iK' 


u + 2K + 2iK' 


snu* = sat 

dnu 


i 

l snu 


—snu 


snu 


—snu 


cnu* = -f'sm 

dnu 


i dnu 

i snu 


— cnu 


—cnu 


cnu 


dnu* = /' ^- 

dnu 


_jcnu 
snu 


dnu 


—dnu 


—dnu 



Table 2 
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We list here all the indentities satisfied by the Jacobian elliptic functions used in 
the course of our calculations 



1 - 


P 


sn 2 usn 2 f 


cn 2 u 




sn 2 u dn 2 v 


1 - 


P 


sn 2 u sn 2 u 


dn 2 v 




;2 2 2 

rsn ucn v 



snu cn« sn(2u) 



(A.5) 



cn 2 u + sn 2 u = 1 ; 
dn 2 u + / 2 sn 2 u = 1 ; 

, . snu cnv dnv ± snv cnu dnu 

sn(u±u) = — ; 

1 — P sn 2 u sn 2 v 

, . cnu cnv =|= snu snu dnu dni> 

cn(u±v) = — ; A.6 

1 — P smusmu 

dnu dnv =F P snu snv cnu cnv 

dn{u±v) = — ; 

1 — P smu smv 

sn(u + v) sn(u - v) = - — — (A.7) 



. . . . U — Oil U. U.11 U . 

cn(u + v) cn(u - v) = — — (A. 



dn(u + v)dn(u - v) = — — — (A.9) 

1 — Psn 2 usn 2 v 



(A.10) 



dnu 1 + dn(2u) 

2 cn 2 u 2 / . v 

sn^u + — 7y- = — — r (A. 11) 

dn 2 u 1 + dn(2u) V ; 

When the module / goes to zero, we have the approximations 
snu ~ sin u — - P{u — sin u cos u) cos u ; 

cnu ~ cos u + -P (u — sin u cos u) sin u ; (A. 12) 

dnu ~ 1 — - P sin 2 u . 
2 

Finally, we also report some useful identities concerning infinite sums of trigono- 
metric functions. These identities are useful for calculations in Sections 3 and 5. 

y sin(2n zDE = i t (A.i3) 

^ 2n-l 4 ' K ; 

n=l 
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n=l 



sin(nx) n — x 



n 



As a consequence of the previous identity we also have 



\n-l 



(-ir * . 

> sm nx 



n=l 



Other useful identities are 



1 11 

V q 2n ~ l cos (2n - l)x = - In - 

^ 2n - 1 v ; 4 1 

n=l 



1 1 + 2q cos x + q 2 



2qcosx + g 2 



— cos(rar) = — - In (l — 2gcosa; + q 2 ) , 



71=1 



r(x) 



T(x + iy) 



k=0 



n i+ 



(x + A;) 2 
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Appendix B 



In this appendix we present the calculations which leads to the exponential repre- 
sentation ( |2.14j ). 

Due to the periodicity of the S-matrix for real values of the rapidity, we can 
always look for its Fourier series expansion with period T of the form 



S(P, a) = — exp 



OO 

2ii\ 



a n e~ nP 



(b.i; 



In order to determine the coefficients a n , let us equating the logarithmic derivative 
of the above expression with the one coming from the S-matrix 



V anne^ = -—E{fi) , (B.2) 

z — <■ n 



where 



sn 



E(B) = sn (2Ko) -7 — t^^t , 7 ' r (B.3) 

K ! 1 ' (sn(«)+sn(2Ka)) (sn («) - sn (2Ka)) 1 J 

and sn'(x) denotes the derivative of the function sn(x). The coefficients a n forn 7^ 

are then given by 

9TC' fl" 

na n = / dt e- 2i $ nt E{t) (B.4) 

Let us consider the complex loop integral 

dte~^ nt E(t) (B.5) 

'c* 

on a closed path C which runs along the boundary of half the fundamental domain, 
i.e. — f < Ret < -|, < Imt < in. Due to the real periodicity of the function 
E(t) on the complex t-plane, it is easy to see that the contributions from the 
vertical lines of the path cancel each others, while the integral on the Imt = in 
line can be rewritten as an integral on the real axis by exploiting the property 
E(t + in) = —E(t). Therefore, for the coefficients a n we have 

na n (l + e 2n ^) = / dt e~^ nt E(t) (B.6) 



nT 
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c 



Inside the domain of integration the function E(t) has two simple poles at t 
ina, Z7r(l — a). Computing the corresponding residues we obtain 



na r 



cosh [^(1- 2a. 



(B.7) 



cosh (zf) 

Inserting this result in eq. ( |B.1| ) and exploiting the parity properties of the expression 
under n — > —n, we can write 



S(P, a) = — exp 



oo 1 

^ n 



cosh ' ra7r2(1 ~ 2a) 



n=l 



cosh (vf) 



sin 



2nn 



p ) + a 



(B.8) 



where ao is the contribution of the zero-mode, which is given by 

nir 2 (l-2a) ' 



cosh 



a = lim 



sm 



(2n^(3\ 2m (5 



T 



T 



(B.9) 



«-Q n cosh (^) 

The zero mode ao ensures in particular that in the limit a^Owe have S(ft, a) — > 1. 
Now, using the identity (|A.15|) to rewrite the zero mode as a series expansion, we 
finally obtain the expression 



S(ft) = — exp 



2zY- 
^ n 



cosh 



( nn 2 (l- 
\ T 



2a) 



■l) n sin 



n=l 



2nix 



ft 



(B.10) 



cosh (ap) 

given in the text, eq. ( |2.14| ). 

In order to obtain the second expression (|2.16|) for the S'-matrix, it is initially 
convenient to split the zero-mode contribution into two series on the odd and even 
numbers, i.e. 



+00/^-1 f2nirft\ 
2i - - sin ' 1 



n=l 



7? 



1 / 

2i > sin 

^ 2n - 1 V 

n=l v 



\ T J 

2{2n - l)n(3 



T 



- 2i 



^ 1 Anvr/A 

n=l v y 



By using now the representation 



— 1 = exp 



n=l v 



27r(2n - 1)/? 
T 
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and adding the term in the exponential of this expression to the zero mode, as a 
result we have that the series of the zero mode on the odd numbers changes its sign 
in front so that we end up with the expression 



cxp 



+00 * 

-2i > — sin 
^ n 



n=l 



2mrp 
T 



(B.11) 



By taking into account the previous Fourier expansion we have then 



exp 



exp 



+00 ^ 

2% V- 
^ n 



cosh 



mr 2 (l-2a) 



n=l 



COSh (2f ) 



+°° 1 sinh s^- sinh . 



— 1 sin 



n(l—a)-K 2 



2mr 
~T~ 



ft 



(B.12) 



n=l 



COSh (?f) 



2mr 
sin I — /3 



i.e. the expression ( |2.16 ) in the text. 
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Appendix C 

In this appendix we report the detailed calculation for the ESG minimal Form 
Factor as given in eq. ( |5.2| ). Let us consider the minimal solution to the monodromy 
equations for the ESG model as given in eq. (|5.1| ). By exploiting the arbitrariness 
of an overall constant in front of F min we can replace —2 sin 2 \nj,pj in eq. (|5.1|) 

( 2rvn 
\ T 



with cos ( "^TjT-f3 ) and take as starting point the following expression 



-Fmin(/?) = exp 



^ ^inhf^isinh'^^ 2 



n=l 



n 



cosh (^f) sinh (^f-) 



2nn ~ 
cos | —=-p 



(c.i; 



Let us first concentrate on the n coefficient in the previous expansion 



sinh ( ) sinh 



n(l — a)7T 2 



(C.2) 



cosh (zf) sinh (^) 
If we use the exponential representation for hyperbolic functions and the series 
expansions 

oo 

(coshx)- 1 = -2 



^ e -(2fc-l)x 



k=l 

oo 



(sinh a;)- 1 = 2 ^ e ~^-i)x 

k=l 



it can be rewritten as 



Cn = ^^(_ 1 )fc+l e -2n[(2fc-l)f 5 ,+(2p-l)^] 
k=l p=l 

x , e -n^(l-2o) _j_ e n^(l-2tt) _ e n^r _ g -n^ 



(C.3) 
(C.4) 

(C.5) 



Now, splitting the fc-sum into a sum on even and odd integers and using the general 
identity 

oo oo oo 

J2J2f( k+ p) = £(* -*)/(*) > 

k=l p=l k=2 

we end up with 

c n = - J> - 1) + + - <£(*) - qUk) - ?{?(*)] , (C.6) 

k=2 
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where we have defined 



-2(2fc-l)£ . n J h \ _ p -2(2fc-2-a)^ 



= e 2(2fc X) r ; q 2 (k) = e 

q3 (k) = e -2(2fc-3+a)^ . ?4 ( fc) = e -2(2fc-l-a)£ . (QJ) 

g 5 (jfc) = e -2(2fc-2+a)^ . ^ = £-2(2^-3)^ _ 

By plugging the expression (|C.6|) of c n into eq. (|C.l|) and performing explicitly the 
sum on n by means of the identity ( |A.17|) , we arrive to the final expression reported 



in the text (eq. ( |5.2| )) 

oo 



fc=2 



where 



W(/3,A;) 



1 - 2q 6 cos (2f/?) + ql) (l - 2q 5 cos (2f/?) + g 5 2 ) (l - 2q 4 cos (2f/?) + gf 
1 - 2 qi cos (2f fi\ + g 2 ) (l - 2q 2 cos (2-0) + g 2 ) (l - 2g 3 cos + qfj 



It is now easy to determine the location of zeros and poles of this function. They 
are simply given by the solutions of the generic equation 1 — 2gcosa; + q 2 = 0, for 
the numerator and denominator factors, respectively. Rewriting this equation as 

(1 - ge*») (1 - qe 4 *) = , (C.9) 

the general solution has the form 

x = ±i\nq + 2mr , neZ . (C.10) 

Applying this result to the factors in eq. (|C.8|) and taking in account the explicit 
expressions for the gj's, we obtain the sets of poles and zeros (|5.4j) , ( |5.5| ) reported 
in the text. 
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Figure Captions 

Figure l.a . Tiling of the /3-plane in terms of the periodic domains of the S- 
matrix. 

Figure l.b . Twists of the elastic branch-cuts of the S-matrix in the s-plane. 

Figure 2 . Analytic structure of the elliptic S'-matrix in its fundamental domain. 
The circles are the zeros of the function and the black circles its poles. 

Figure 3 . Phase-shift (in units of n) of the Elliptic Sinh-Gordon S'-matrix for 
particular values of the parameters (/ = 0.25 and a = 0.3). 

Figure 4 . Phase-shift (in units of n) of the Elliptic Ising Model with the period 
relative to I = 0.25. 

Figure 5 . Poles and zeros in the /3-plane for a = S — i^. They are split by 25 
and annihilate each other for 5^0. In this figure the rectangle delimitated 
by dashed lines is taken as the fundamental domain. 

Figure 6 . Diagram associated to the pole singularity of the S-matrix near the 
resonance R n . 

Figure 7 . Plots of the two-point correlation function of the field Q(x) versus 
M.R. The full line is relative to the usual thermal Ising model (/ = 0), the 
short-dashed line to the value I = 0.25 and the long-dashed line to I = 0.75. 

Figure 8 . Central charge c versus the modulus / for the Elliptic Ising model 
obtained by the c-theorem sum-rule. 

Figure 9 . Anomalous dimension of a(x) and n(x) versus the modulus / for the 
Elliptic Ising model, obtained by the 77 sum rule. 
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Figure 10 . Analytic structure of F min ((3) of the ESG model along the imaginary 
direction placed at f3 — nT . The circles are the zeros of the function whereas 
the black circles represent its poles. 

Figure 11 . Graphical multiplication of F min (/3) by F m - m (f3 + in). 
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